n = 5;

b = 25.867325;

A = [ ... 
  1.075334,  0.526197,  -3.608734,  0.657207,  0.990262,
  0.526197,  -0.867184,  3.377548,  3.454253,  1.561950,
  -3.608734,  3.377548,  1.450808,  1.426643,  1.431982,
  0.657207,  3.454253,  1.426643,  2.818069,  3.047428,
  0.990262,  1.561950,  1.431982,  3.047428,  0.977788,
];

C = eye(n);
